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Abstract 

We present observations of the J=l-0, 2-1, and 3-2 rotational transitions of N2H+ and N2D+ 
towards a sample of prototypical dark clouds. The data have been interpreted using non-local radiative 
transfer models. For all sources previously studied through millimeter continuum observations, we find 
a good agreement between the volume density estimated from our N2H"'" data and that estimated from 
' the dust emission. This confirms that N2H"'' depletion is not very efficient in dark clouds for densities 

' as large as 10^ cm~^, and also points out that a simultaneous analysis based on mm continuum, N2H"'" 

and N2D"'" observations should lead to rehable estimates for the temperature and density structure of 
cold dark clouds. 

^ I From multi-line modeling of N2H'^ and N2D"'", we derive the deuterium enrichment in the observed 

^ . clouds. Our estimates are similar or higher than previous ones. The diflFerences can be explained by 

the assumptions made on the cloud density profile, and by the chemical fractionation occurring in the 
I clouds. For two of the observed objects, L183 and TMC2, multi-position observations have allowed 

us to derive the variation of the N2D"'"/N2H"'" abundance ratio with the radius. We have found that it 
decreases by an order of magnitude for radii greater than a few 0.01 pc (i.e. outside the central cores). 
Inside the dense condensations, the fractionation is efficient and, compared to the abundance ratio 
expected from statistical considerations based on the cosmic D/H ratio, the deuterium enrichment is 
O ■ estimated to be ~ 0.1-0.5 10^. 

An important result from our observations and models is that the interpretation of deuterated 



o 



^ . molecular species emerging intensities in cold dark clouds requires specific radiative transfer modehng 

^ ' because the excitation conditions for the deuterated species can be quite diff'erent from those of the 

main isotopologue. Moreover, the use of three rotational fines for N2H+ and N2D+ allows to constrain 
the size of the emitting regions for each species and to determine accurately the volume density. This 
J> ' enables to draw a detailed picture of the spatial variation of deuterium enrichment. 

CN , Subject headings: line: formation : profiles — molecular processes — radiative transfer — ISM: clouds 

O . : molecules : abundances — ISM: individual (TMCl, L183) 

p , 

\^ . 1. INTRODUCTION Moreover, from NH3 inversion lines it has been found 

O ; Over the past decades, important efforts have been de- that the clouds are nearly isothermal | |Benson fc Mye rs| 

, voted to improve our understanding of the physical con- |198ai Tafalla et al. || 2004D . Such observational constraints 

O . ditions in protostellar clouds. These studies have en- ^re of great interest to understand the physical mecha- 

abled to draw an evolutionary sequence of the stages "^^^^^ t^^^t '^"^e the collapse. However, the balance be- 

prior to the formation of protostars, through the mor- ^ween gravitational energy and support mechanisms such 

phological, dynamical and chemical characteristics of the ^s thermal p ressure, turbu ence o r magnetism is yet to be 

^ clouds. Recent progress in detection capabilities in the understood ( |Aikawa et al. || 2005| ). Moreover, theoretical 

■•^ . far infrared, sub-mm and mm windows, have permitted '^t^'^ie'^ predict that the initial conditions strongly affect 

to constrain the density distribution of molecular clouds the way the collapse of the cloud evo ves. In particular 

through measurements of dust absorption and emission dependent studies, which couple both dynamical 

(see e.g. lWard-Thomnson et al.l[T994l . [T999h . It has been f^^ chemical processes, have shown that molecules can 

found that the geometrical structure of cold dark clouds ^e used as tools to probe the stage of evolution prior 

is consistent with an inner region of nearly constant den- ^^"^ formation of protostars. These studies also pre- 

sity and a surrounding envelope that has a density which ^ict that deuterium fractionation increases with time and 

decreases as a power law. The density distribution thus that the D /H ratio can reac h high values for a variety of 

obtained is often used as a starting point for molecu- species ( |Roberts et al. || 2003| j. 

lar fine studies and allows to const rain the spat ial dis- The goal of this work is to put observational constraints 

tribution of the molecules (see e.g. iTafaliret all[200l . ^^e spatial structure of a sample of dark clouds 

through observations of a late type molecule, N2H'^, and 
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sential informations on molecular physics data, such as 
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collisional rate coefficients, line intensities and frequen- 
cies. Molecules often used as tools, like N2H+, have been 
interpreted so f ar with crude estimates of these param- 
eters. Recently, iDaniel et al.l (|2005f ) have computed col- 
Hsional rate coefficients between N2H+ and He for the 
range of temperatures prevaili ng in cold dark clou ds and 
low mass star formi ng regioii s. iDaniel et al.l (|2006f ) (here- 
after referred to as iPaper J ) have modeled, using these 
new rates, the emission of N2H"'" in dark clouds. A wide 
range of physical conditions has been explored in order to 
have a general view of the different processes that could 
lead to the observed intensities of N2H+ hyperfine lines. 

Observations are described in section [21 In section [3] 
we present the results obtained for each cloud and we 
derive the physical conditions of the gas as well as N2H''" 
and N2D+ abundances. The conclusions are given in sec- 
tion [H The excitation processes occurring for deuterated 
species are analyzed in the Appendix. 

2. OBSERVATIONS 

N2H"'" and N2D'^ observations were performed during 
several observing sessions. Addit ionally, we make use of 
N2D+ data already discussed by lTine et al.l 1)20001 ). 

Observations of the N2D+ J=2-l Hne have been ob- 
tained at the IRAM-30m telescope (Granada, Spain) in 
July 2000 during average weather conditions. The point- 
ing was checked regularly on nearby continuum sources, 
and the focus was checked on small planets. The data 
have been taken in frequency switching mode, with the 
mixers connected to the versatile, high spectral resolu- 
tion, correlator VESPA. The N2D+ J=l-0 data have 
been taken in August 2004 when the tuning range of the 
3mm mixers was extended to frequencies lower than 80 
GHz. Simultaneous observations of N2D"'" J=3-2 were 
performed using two mixers tuned at this frequency. The 
system temperature was ~ 250 K at 77 GHz and ~ 500 
K at 231 GHz. Particular care was given to the inten- 
sity calibration as the mixer tuning is done without side 
band rejection at these low frequencies. The image side 
band gain has been measured on an offset position or a 
continuum source for each object, and applied during the 
data processing. 

Further observations of the J=l-0, 2-1 and 3-2 lines 
of N2H+ and of the N2D+ J=3-2 Hne were carried out 
in January 2005 with the IRAM-30m telescope. Four 
receivers tuned at the frequency of these lines were used 
simultaneously. Pointing was checked every hour on 
continuum sources close to the targets. The relative 
alignment of the receivers was checked with Saturn and 
Jupiter and found to be better than 2.0". Our observa- 
tions indicate that the offsets on the sky from the point- 
ing derived with the 3 mm receiver were 1.0", 2.0" and 
1.0" for the 2 mm, 1 mm low frequency and 1 mm high 
frequency receivers respectively. Focus was checked at 
the beginning and end of each observing shift. It was 
found to be stable for the entire observing run. 

The weather was excellent during observations with a 
typical water vapor amount above the telescope site of ~ 
1.0 mm, or even better. Calibration was performed using 
two absorbers at ambient and li quid nitroge n tempera- 
tures . The ATM package l|Cernic haro 1985; Pardo etldl 
I2OOII I was used to estimate the significant atmospheric 
parameters for the observations. The typical system tem- 
peratures were 100 K, 1000 K and 250 K at 3, 2 and 1 mm 



respectively. The J=2-l line was observed in upper side 
band which avoided the possibility of rejecting the image 
side band in the 30-m receivers. The large system tem- 
perature for this transition is due to the high opacity of 
the atmosphere at 186 GHz and to the fact that the fre- 
quency is at the edge of the tuning range of the receiver. 
Calibration uncertainties for the J=l-0 and 3-2 lines are 
well below 10% (the standard calibration accuracy under 
good weather conditions). However, the J=2-l Hne has a 
larger uncertainty due the effects commented above, and, 
in order to have an estimate of the calibration accuracy 
we observed the continuum emission of Jupiter, Saturn 
and Mars at different elevations. We have found that the 
intensity of the J=2-l transition was always compatible 
with the expected one. However, these double side band 
observations of the planets are not very sensitive to the 
side band rejection. Nevertheless, these data allow us to 
derive beam efficiencies for each frequency. They are in 
excellent agreement with those published for the 30-m 
IRAM telescope. For some strong sources we observed 
several times a reference position and we noted that the 
J=2-l line intensity had calibration errors up to 20% de- 
pending on the elevation of the source. We conclude that 
with the configuration adopted for the 2 mm receiver, the 
upper signal side band has a gain of 0.35 rather than 0.5 
(perfect DSB). This fact introduces systematic errors in 
the determination of the sky opacities. It is impossi- 
ble to estimate the actual error for each source. Hence, 
for the J=2-l line we assume a typical caHbration error 
of 20%. Nevertheless, most observed J=2-l intensities 
agree with those expected from the results of the J=l- 
and 3-2 lines (see Sect. Only for TMC2, L1489 
and L1251C it has been necessary to correct the J=2-l 
intensities by multiplicative factors of, respectively, 0.8, 
0.8 and 0.85 in order to have data consistent with the 
well calibrated J=l-0 and 3-2 fines. 

The observations were performed in frequency switch- 
ing mode. As all sources have narrow lines, the baseline 
resulting from folding the data has been removed by fit- 
ting a polynomial of degree 2. Figure [T] shows the N2H+ 
and N2D+ observed lines in selected sources and Table 
[1] gives the position of the clouds observed in this paper 
(see section [S]). The intensity scale is antenna tempera- 
ture, T^. In order to compare the results of our models 
with observations, we have to convolve model results with 
the beam pattern of the 30-m IRAM telescope. We have 
adopted half power beam widths (HPBW) of 27", 13.5" 
and 9", beam efficiencies of 0.76, 0.59 and 0.42, and error 
beams of 350", 220", and 160" for the J=l-0, 2-1 and 
3-2 lines of N2H+. For N2D+, we adopted HPBW's of 
32", 17" and 10.5", beam efficiencies of 0.79, 0.69 and 0.52 
and the same error beams as for N2H"'' lines. The whole 
data set for all observed sources is shown in Figures 3 to 
15. 

3. RESULTS 

Over the past decades, mm continuum observations 
and star count analysis (|Ward-Thompson et al.l Il994[ ) 
have revealed that the density structure of cold dense 
cores closely matches a Bonnor-Ebert sphere, with an 
inner region of nearly uniform density and steep outer 
edges where the density decreases as t~p with p = 2 - 
2.5. This description is valid for r < 0.1 pc typically, 
and at larger scales (0.1 - 1 pc), the density structure of 
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the outer envelope has be en shown to vary f rom s ources 
to sources. For example, ICernicharo et al.1 (|l985f ) have 
found that a shallow density profile, n(r) oc r~^/^, is ad- 
equate to model the density structure of the envelope of 
the Taurus-Auriga complex. In this section, we use a 
density profile which corresponds to n(r) = no for r < ro 
and n(r) = no(ro/r)^ for r > rp. In order to characterize 
the clouds, we determine the value of the central density, 
i.e. no, and the radius of the inner fiat region, i.e. ro. For 
L1517B, we compare the results obtained using this den- 
sit y profile with the p henomenological density law used 
in lTafalla etHI (|2004l l. 

The observed line profiles for N2H"'' and N2D+ J=l- 
0, 2-1, and 3-2 fines are shown in Figures 3 to 15. In 
order to model our data we have found that some ad- 
ditional constraints were necessary. In particular the 
extent of the emission has been found to be a criti- 
cal parameter to derive physical conditions. Our data 
are in most cases related to one single observed posi- 
tion (except for L183 for which we have a map in sev- 
eral transitions). Fortunately, some clouds presented 
in this section have been mapped in th e J=l-0 tran- 
sitions of N2H+ by ICa selli et iaU l|2002l l. iTafalla et all 
l|2004| l. and lCrapsi et al . (2005) (In the discussion of the 
sources, if not spec ified, N2H+ maps refer to data from 
ICaselli et all (|2002h ). For six sources (L63, L43, L1489, 
L1498, L1517B, and TMC-2) they have reported J=l-0 
integrated intensity versus impact parameter. Thus, we 
have checked that the profile of the J=l-0 integrated in- 
tensity provided by our models was consistent with those 
data. For the other clouds, similar data could not be 
found in the literature. 

The radiative transfer models aim at reproducing our 
J=l-0, 2-1 and 3-2 line intensities and the J=l-0 ra- 
dial emission profile. A discussion concerning the as- 
sumptions made to model N2D+ is given in the Appendix 
(that deals in particular with the determination of col- 
lisional rate coefficients for the isotopologue of N2H"'"). 
Observed positions, central line intensities and adopted 
distances to the sources are given in Table [TJ We re- 
port in the same table the J=l-0 emission half power 
radius (noted r). This corresponds to the radius which 
encloses half of the J=l-0 total intensity which emerges 
of the cloud. The convolution of the emerging line in- 
tensities with the telescope beam pattern has been done 
taking into account the characteristics of 30-m IRAM 
radio telescope (see section 2) and those of the 14-m 
FCRAO antenna. For the latter we have used a HPBW 
of 54" at 91 GHz and a beam efficiency of 0.51. Never- 
theless, the error beam size for the FCRAO antenna is 
not known but must be taken into account. In our mod- 
eling, we arbitrarily assume a size of 180". As quoted 
above, the radial intensity profile enables to constrain 
the way density and abundance vary through the cloud. 
We stress that the density profiles determined in the fol- 
lowing section depend essentially on the observations of 
the J=l-0 line. Hence, an accurate description of the 
radiation pattern of the FCRAO antenna is needed to 
obtain the best estimates for the physical parameters of 
the clouds. We have checked the effect of introducing 
different error beam sizes for the FCRAO antenna, and 
the largest differences are obtained when the error beam 
is not taken into account, which leads to J=l-0 line in- 



tensities 20-30% lower for the same physical conditions 
and source structure. This effect could be compensated 
by increasing the size of the central core, ro, by 5"-10", 
decreasing its density, no, and increasing the abundance 
of N2H+ in the outer region of the clouds. 

For some sources, we have included infall velocity fields 
in order to reproduce the spectral features seen in the 
J=2-l line. The profiles used consist of a w oc r func- 
tion which corresponds to infalling cores for all sources, 
except for L1517B. Despite the inclusion of these veloc- 
ity fields, we have observed frequency offsets between the 
observations and the models. In table HI we report the 
"Vlsr determined from our J=l-0 data, the corrections 
to this value for the other lines (Vlsr and corrections 
are determined by a x square fitting), and the velocity 
gradient found in the sources. For all of them, we find 
that the J=2-l and J=3-2 lines have \lsr ~ 30 m s~^ 
larger than the \lsr of the J=l-0 line. This would be 
consistent with the view of a static core with an infalling 
envelope as the J=2-l and J=3-2 fines arise from the 
densest part of the cloud. However, we think that the 
effect, 20 and 30 kHz for the J=2-l and J=3-2 fines 
respectively, is within the a ccuracy of the labor atory fre- 
quencies of N2H+ (see e.g. [A mano et al.ll2005l l. Conse- 
quently, the range of model parameters is limited by the 
laboratory accuracy in the frequencies. 

3.1. Individual core properties 

In this section, the cloud physical parameters (see 
Tables [2] and [3]), obtained from the modefing of our 
data, are discussed. In the following discussion, we 
use the notation ro = X(N2D+)/X(N2H+) to refer to 
the isotopologues abundance ratio. This quantity dif- 
fers from the column density ratio (referred to as i? = 
N(N2D+)/N(N2H+)) in the case of three clouds (L63, 
TMC2 and the two molecular peaks in L183) for which 
rj) has been found to change with radius. R represents 
in all cases an averaged value along the line of sight of 
the N2H"'" volume density. 

For each cloud, the error bars quoted for the physical 
parameters are derived by use of a reduced ■ These are 
obtained by assuming that the parameters of the models 
are independent and thus, it does not account for the 
correlation among them. Nevertheless, it gives an indi- 
cation on the sensitivity of the emerging spectra to the 
model's parameters. The same criterion was used for 
both N2H'^ and N2D+ spectra. The uncertainty on ro 
is then found important due to the lower signal to noise 
ratio of the N2D"'" data and this leads to large error bars 
for the N2D+ column density (varying from 20% to 100% 
depending on the source) while typically, the N2H+ col- 
umn density is estimated with an accuracy of 10-20 %. 

3.1.1. L1517B 

Modeled fines and observations of the J=l-0 and J=2- 
1 transitions of N2H''" and N2D"'' are shown on figure 
[3l When comparing the line profiles of the two isotopo- 
logues, it is important to note the differences which are 
introduc ed by th eir respective abundances, i.e. opaci- 
ties (see iPaper J ) . We see that the intensity ratio be- 
tween the different J=l-0 hyperfine components are best 
reproduced for N2D"'" which is optically thinner than 
N2H"'". This fact supports the idea that hyperfine inten- 
sity anomalies are due to radiative effects which occur at 
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high opacities. For the J=2-l Hnes of the two isotopo- 
logues, we see that the main central emission component 
is self absorbed for N2 H"'" while it is not for N2D"'". More- 
over, as commented in lPaper"3 . the line profile seen in the 
J=2-l line of N2H"'" is characteristic of expanding flows. 

A detailed study of L1517B has been reported in 
iTafalla et al.l (|2004f ) where the structure of the cloud is 
derived from maps of the NH3 inversion line and 1.2 mm 
continuum emission. The position we ob served corre- 
sponds to the N2H"'" peak intensity given bv lCaselli et al.l 
|200a), which is situated at (+15",+15") of the N2H+ 
emission pea k in th e higher resolution map reported by 
iTafalla et al.l (|2Q04f ). The latter position is also a max- 
imum of the 1.2 mm emission map. In our models we 
to ok into account the J =1-0 and J=3-2 lines reported 
by lTafalla et al.1 l|2004f ) towards the continuum peak. On 
figure O we compare the J =1-0 radial integrated inten- 
sity of our model to the observational data reported in 
the latter study. 

To characterize the velocity structure of L1517B, we in- 
troduced a linear outward gradient o f 1.4 km s~^ pc^^ for 
radii beyond 20". As pointed out in lPaper"3 . it is rather 
difficult to constrain velocity fields from observations of 
the J=l-0 and J=3-2 lines of N2H'^, while the larger 
opacity of the J=2-l line enables a better estimate. 
rVom observations of CS emission Hnes, JCafalla et aQ 
l|2Q04l l derived a value for the outward gradient of 2 km 
s~^ pc~^, for radii larger than ~ 70". Hence, this velocity 
field is consistent with their observations of the J=l-0 
and J =3-2 Hnes of N2H+ observed towards the central 
position. 

On Figure [21 we compare the J=l-0 radial intensities 
obtained with two types of density profiles. The first is 
the one quoted above (labeled (1) in Table [2|) and the 
second corresponds to the phenomenological profile used 
bv .Tafalla et al. (200J) (labeled as (2)). Using the profile 
(2), we obtain the same estimate than in the latter study, 
i.e. no = 2.0 10^ cm ~^ and ro = 35". Moreover, we see 
that with the present observational data, we are rather 
insensitive to the exact form of the function which de- 
scribes the density variations in the cloud, since the two 
density profiles give similar results for both line shapes 
and the J=l-0 integrated intensity profile. 

Using the two density profiles mentioned above, we also 
derive different estimates for the N2H"'" abundance, i.e. 
1.6 10-1° with (1) and 2.5 10-^° with (2). In the latter 
case, our estimate is higher by ^ 70% with respect to 
the one obtained bv ITafalla etaH (|2004[ 1. i.e. X(N2H+) 
= 1.5 10"^°. The difference is due to the use of different 
sets of coHisional rate coefficients in each studies. 

We have determined the D/H ratio that allows to re- 
produce the J=2-l line of N2D+ for the two density pro- 
files. For the density profile (1), we derive R = 0.12to o3 
which is a facto r 2 higher than the estimate of 0.06 ± 
0.01 obtained bv lCrapsi et al.l (|2005f ) toward the density 
p eak position. 

iBenson fc Myers! (|l989h have estimated a mass of 0.33 
Mq within a radius of 1.4', which corresponds to a mean 
volume density of n(H2)=6 10^ cm"^. For the same ra- 
dius, we obtain a mass of 2.5 M©. On the other hand, 
the analysis of e mission maps at 850 fim and 450 fim 
ijKirk et al.ll2005f ) lead to an estimate of n(H2) ~ 4-5 10^ 
cm~^, in qualitative agreement with the central density 
assessed from our N2H+ data. 



3.1.2. L183 

L183 is a well studied object due to its proximity 
(d=110 pc). It has a non spherical shape with a north- 
south axis in the density distribution. Recent stud- 
ies point out a complex chemistry in the cloud, with 
four main molecular peaks showing characteristics which 
are interpreted as signs of different stages of evolution 
ijPickens et al.i r2000l. The N2H+ maximum of intensity 
is situated near the densest part of the cloud (molecu- 
lar peak C). A secondary maximum situated ~ 3' to the 
north (molecular peak N) is also seen in the N2H+ in- 
tensity map. For this source, the observed position for 
molecula r peak C cor r espon ds to the intensity peak as 
given bv ICaseHi etldl ||2002D . It is placed at (-5",+15") 
with respect to the peak positi o n of th e hig her resolution 
maps given by iDickens et al.l l|2000l ) and iPagani et al.1 

(Hop). 

Figure [4] shows peak C observations and models of the 
J=l-0, 2-1 and 3-2 transitions of N2H+, and of the 
J=3-2 transition of N2D''". Figure [5] shows the J=l- 
and 3-2 Hnes of N2H+ and the J=l-0, 2-1 and 3-2 
Hnes of N2D+ for peak N. Figures [6] and [7] present maps 
of the J=l-0 and 2-1 Hnes of N2D+ obtained towards 
pe ak C, the latter b eing obtained with the same data as 
in lTine et all l|2000i ). 

There are a number of stu dies regard i ng the dust emis- 
sion pro perties of L18 3 (e.g. iKir k et al] ||2005D at 450 and 
850 um. lJuvela et al. (2002), Pa gani et alJ (|2003f ) at 200 
/zm and IPag ani et al.l (|2005ll at 1.2 mm). From these 
studies, it seems that the region of molecular peak C is 
associated with two continuum sources FIRl (south) and 
FIR2 (north) separated by - 1.5'. The region of FIRl 
and FIR2 corresponds to the maxir num extinction an d 
minimum temperature in the cloud ijJuvela et al.ll200"^ . 
Also, there might be a gradient of temperature in this 
region or changes in the composition of the grains since 
FIRl is detected at 200 and 450 ^m but remains unde- 
tected at 850 /im. Moreover, the emission peak FIR2 at 
850 fim is shi fted to the North compared to that at 200 
and 450 ^m IjKirk et al.ll2005[ l. while the 1.2 mm peak 
position is shifted to the South. It is worth noting that 
the maxima of the 1.2 mm and N2H"'" maps reported in 
IPagani et al.l (|2005f ) do not coincide. From this study, it 
appears that the maxima of these two maps are shifted 
by ~ 25" and that the H2 volume density reaches 2 10^ 
cm~^ at the position of the 1.2 mm map emission peak. 
For simplicity, we have modeled the observations of peak 
C assuming that the density center of our model is at 
the position of the N2H"'' emission peak (i.e. offset by 
(+5", -15") with respect to the coordinates given in table 

m- 

For the spectra belonging to peak N, our observations 
are off b y (+2",-40") with resp ect to the reference position 
given by lDickens et al.l (|2000h . We have assumed in our 
models that the density peak is centered at the position 
of our observations which corresponds to the secondary 
maximum in the N2H"*" map. 

For peak C, in order to fit the observed asymmetry in 
the J=2-l transition, we have included a veloci ty gradi- 
ent of Av = -0.2 km s"^ pc"^. As discussed in lPaperl . 
this Hue is sensitive to small velocity gradients, due to 
the presence of multiple hyperfine components in a small 
frequency interval. For peak N, no velocity gradient was 
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included since, for this position, the observed lines are 
insensitive to such small velocity fields (the J=2-l line 
was not observed). It is worth noting that the J =2-1 
transition of N2H"'" obtained towards peak C has a pecu- 
liar shape that is not observed in any other cloud of the 
present sample. Hence, the abundance profile obtained 
for this source is unique, among all our models, since we 
find that the best agreement between models and obser- 
vations is obtained by increasing the abundance at large 
radii. 

For the two peaks, we obtain similar N2H~'' abun- 
dances, but different density profiles: peak C is more 
dense and centrally peaked than peak N. For the former 
position, our density estimate is consistent with that ob- 
tained from sub-m m observations, i.e. n(H2) = 10^ cm~'^ 
l|Kirk et al.ll2005l l. 

In order to model the N2D"'" map shown in Figures 
[6] and [3 the relative abundance has been set to rn = 
0.43 for radii below r ~ 0.021 pc. For larger radii, the 
relative abundance is lowered to 0.05. Note that this de- 
crease in abundance is necessary to reproduce the line 
shapes as it reduces the self-absorption induced by the 
external layers of the cloud, in the J=2-l main compo- 
nent of N2D'^. On the contrary, the main component of 
the J=2-l transition of N2H"'" is strongly self-absorbed 
and is better reproduced by increasing the N2H+ abun- 
dance at the outermost radii. Since our model considers 
spherical geometry while the morphology of the cloud is 
extended along a North-South axis, we have failed to re- 
produce the J=2-l transition for the whole map. Hence, 
the intensities are well reproduced for declinations be- 
tween S = -20" and +20" but are underestimated along 
the North-South axis. Toward the position given in ta- 
ble [H we find a column density ra tio R = 0.33ln;|i w hich 
is similar to the one obtained by iTine et all (|2000h at a 
position offset by 5" to the north, i.e. R = 0.35. This 
value is higher than the one obtained by ICrapsi et al.1 
l|2005D at a position offset by (-5",-15"), i.e. R = 0.22 
±0.04. 

From the observations towards peak N, we obtain val- 
ues for the deuterium enrichment of r_D = 0.76 for radii 
below 5 10"'^ pc and 0.05 beyond. Note that the radius 
which encloses the central D-rich region is determined 
from a single position analysis. Because the angular size 
of this central region is intermediate between the HPBW 
of the J=l-0 and 3-2 lines of N2D+, (32" and 10" re- 
spectively), the J=3-2 line is favored compared to the 
J=l-0 one. In order to better constrain the variation 
of the D/H ratio with radius, observations at different 
positions are needed. For peak N, we derive a column 
density ratio of R = 0.15 . 

Recently, a similar work bas ed on N2H"'" and N 2D"'" 
observations was carried out bv iPaganT et all l|2007l l for 
molecular peak C. This latter study differs from ours in 
the location of the center of the density profile: it was 
set coincident with the maximum of the 1.2 mm map 
which is shifted 25" to the south compared to that as- 
sumed in the presen t work. Moreover, the ID radiative 
transfer code used in lPagani et al.l (|2007f ) work treats line 
overlap and the convolution with the telescope beam was 
made using an approximation that takes into account the 
north-south extent of the cloud. At the position of the 1.2 
mm map peak, they found that the temperature ranges 



from 7K (from N2H+ data) to 8K (from N2D+ data) and 
that the maximum H2 density is 3 times higher than that 
derived in the present work. Note that the density profile 
derive d here and the "best model" profile in lPagani et al.l 
{^OOf) are similar beyond 30": the latter corresponds to 
densities which are higher by ~ 25%. Th e main differ- 
ence c oncerns the N2H"'' abundance profile. iPagani et al.l 
l|2007l l have found that X(N2H+) drops by a factor 6^3^ 
in the inner 20" region. This decrease is interpreted as 
depletion and it becomes efficient at densites n(H2) = 
5-7 10^ cm~^. At our position, we do not see evidence 
for such a strong decrease in abundance: the abundance 
we find is ~ 1.6 time lower below 40" where the den- 
sity reach n(H2) = 2 10^ cm~'^. Given that the source 
strongly departs from spherical geometry, this variation 
could be a consequence of a mis-interpretation of the 
real excitation conditions and prevents us to interpret 
this result as a sig n of depletion . Con sidering the deu- 
terium enrichment. iPagani et al.1 (|2007l ) find that at the 
1.2 mm map peak, ro varies from 0.7 ±0.12 (inner 6") to 
0.05 (beyond 40"), which is consistent with the present 
findings. 

3.1.3. TMC2 

In the present study, we report N2H+ and N2D+ obser- 
vations for two positions towards TMC2 which are sepa- 
rated by 40". Figure [8] shows the J=l-0, 2-1 and 3-2 
observed lines of N2H"'", and the J=3-2 line of N2D"'' for 
the central position. For the (-33",-23") position, N2H"'" 
J=l-0 and N2D"'' J=2-l fines are shown. Figure[2] shows 
the J=l-0 radial integrated inte nsity of the ni o del co m- 
pared to observational data from ICaselli et al.l (|2002h . 

This source has been mapped in ma ny molecular 
lines ( e.g. C^^O. D CQ+ and H^'^ CO+ by [ButnerelaLl 
IT995I1. H"CO+ by lOnishi et al.1 (1200^ and C^^O by 
lOnishi et al. lll996ll') as well as in 1.2 mm continuum 
emission bv ICrapsi et al.1 l|2005f ). Large scale visual ex- 
tinction maps covering the w hole region are reported by 
ICernicharo et al.] (|l984l . [l985f ). Contrary to other clouds, 
N2H"'" and 1.2 mm maps seem to be uncorrelated: while 
the N2H+ peak has a bright 1.2 mm emission counter- 
part, there are several maxima in the 1.2 mm map which 
do not correspond to any N2H"'" peak. In this work, we 
derive an average central volume density of n(H2) ~ 3.5 
10^ cm~^, which is in good agreement with the value of 3 
10^ cm~'^ derived from the analysis of the 1.2 mm emis- 
sion. However, this value is higher than the estimates 
obtained from other molecular species. The discrep- 
ancy partly arises from a different spatial l ocation of the 
molecules (see the map of C^^O reported by'CraDsi et alj 
Hooa) and the map of Hi3C0+ of Onishi et al."(2002)V 
On t he other hand, N2H+ and NH3 lines (Myers et edl 
I1979I I seem to arise from the same region and an average 
volume density of 2.5 10^ cm~^ was derived from the lat- 
ter mol ecule. Similar values for the average density are 
given in lButner eraP lfT995) . based on C^^O and DCO+ 
lines. Nevertheless, high er density estimate s have been 
obtained using Hi3C0+ (jOnishi et al.ll2002l l. The map 
of this molecule shows that two emission maxima are lo- 
cated within the NH3 map corresponding to regions of 
average volume density ~ 10^ cm~^. 

From NH3 lines, the mass enclosed within a radius of 
3.6' h as been estiraated to be ~ 16 M0, and from C^^O 
lines (jOnishi et al.ll20d^ . the mass enclosed within 7.8' 



6 



is estimated to be ~ 50 M©. The density profile given in 
table [2] corresponds to enclosed masses of 50 Mq and 121 
Mq for the same radii. We stress that the value adopted 
for ro is the most important parameter to derive the mass 
of the cloud and, in order to obtain similar values than 
those quoted above and using the same central density, 
we would have to adopt a value of ro ~ 35". Adopting 
such a small value for ro does not allow to reproduce the 
J=l-0 radial integrated intensity. 

To reproduce the N2D"'" spectra obtained towards the 
two positions, we have decreased the abundance ratio by 
a factor 7 at radii > 0.024 pc. It is worth noting that 
a single position analysis enables to reproduce indepen- 
dently each one of the two N2D"'" spectra, assuming a 
constant value for ro. For the central position, such an 
analysis leads to a ratio of rjj ~ 0.25. 

3.1.4. TMCl(NHz) 

Figure [9] shows observations of N2H+ and N2D+ to- 
wards two positions separated by ~ 35". The reference 
position is given in table [1] and the second position is 
located at (-35", +9"). Towards both positions, two ve- 
locity components are identified in the observations of 
the J=l-0 line of N2H"'". In order to model the faintest 
component, we have simply added the contribution of 
two emitting regions. This way of computing the emerg- 
ing spectra is motivated by the fact that the difference 
in the N lsr of the two components is larger than the 
Hnewidth: the determined \ lsr are 5.98 km and 
5.68 km , the former being the most intense one. 
For both positions, the faintest velocity components have 
been determined to be similar. 

The TMCl ridge is known to have a complex velocity 
structure with a velocity gradient of ^ 0.2 km s~^, per- 
pen dicular to the major axis of the ridge (jOlano et al.l 
119881 , and cores at different Vlsr for the same line 
of sight. Particularly, two components at ~ 5.6 km 
and 6.0 km have been ident ified around the 
cyanopolyyne peak ijTolle et al.l llQSlf ). Also, around 
the ammonia peak, two possible components at Vlsa 
^ 5.9 km and 5.3-5.4 km are reported by 
iHirahara et al.l l|l992l l from observations of C'^^S tran- 
sitions, the latter being identified as a low density core. 
The two Vlsr components identified in this work are 
also seen in SO observations ijLique et al.ll200^ . 

As for the previous sources, we obtain a higher den- 
sity estimate than in previous studies based on the emis- 
sion from other molecules. See for example the detailed 
analysis of the density st r ucture along the TMCl ridge 
reported bv iPratap et al.1 (|1997[ ). At the position of the 
N2H+ peak (ammonia reference position), the density is 
estimated to be 10^ cm~^ from the analysis of several 
HC3N lines. Note that the structure of the cyanopolyyne 
peak has been probed th rough interferometric observa- 
tions l|Langer et al.|[T995l l. It has been found that the 
core of the cloud consists of several high density con- 
densations (n(H2) ~ 3 10^ cm~^) with diameters rang- 
ing from 10" to 30", which are embedded in a lower 
density medium. This result was further confirmed by 
iToth et all (|2004l l who analyzed ISOPHOT images. 

From the J=2-l line of N2D+, we derive R ~ O.lOt^;^^^ 
for the ammonia reference position, which is in good 
agreem ent with the value of 0.08 reported bv lTine et al.l 
l|2000t) . 



3.1.5. LI498 

Figure [10] shows the J=l-0 line of N2H+, and the 
J=l-0 and 2-1 lines of N2D+. L1498 is a nearby 
molecular cloud located in the Tau rus com plex. This 
cloud has been observed in H2CO llYoung et al. 2004), 
CS (Tafalla et al."2004D and N2H+ (jTafalla et al.ii200l : 



iShi rlev et al. 2005) and its physical pr operties have 
been derived from ISOPHOT obser vations (iLanger et al.1 
I2OOII I and millimeter observations ijShirlev et al.ll2005[ ). 

As for L1517B , a det ailed study of L1498 has been done 
bv lTafalla et all (|2004f ). and we have used this work as a 
starting point for our models. The position we observed 
in this source is at (+10",0") of the 1.2 mm emission peak. 
In our model, we have assumed that the center of the den- 
sity profile corresponds to the continuum peak. Figure [2] 
shows the J=l-0 radial integrated i ntensity comp ared to 
the observational data reported bv lTafalla et al.l |[2004). 
As noted in that study, a decrease of the abundance in 
the outer envelope is needed to reproduce the J=l-0 ra- 
dial intensity profile. Otherwise, the intensity would be 
overestimated at large radii. 

The values we have derived for the density and for 
the N2H"'" abun d ance are consistent with those found 
by iTafalla et aP (|2004l l. As for L1517B, we derive a 
higher abundance than in the latter study (X(N2H+) = 
2.5 10-1° compared to 1.7 10"^°). This is related to 
the use of different collisional rate coefficients and den- 
sity profiles in each study. Moreover, we derive the same 
central density, but with a smaller value for ro due to 
the different analytical expressions of the density profile. 
Nevertheless, the central density reported in this work 
is higher than what is usually found for this cloud. For 
example, IShirley et al.l (|2005( i reported a density ~ 1-3 
10"* cm"^ from continuum observations at 350 /im, 850 
/zm and 1.2 mm. 

The value we have found for the column density ratio, 
i.e. R = 0.0 7j^Q 03, is larg er by a factor two than the one 
reported bv lCrapsi et al.l (|2005f ) at the same position. 

3.1.6. L63 

Figure [11] presents the J=l-0, 2-1 and 3-2 lines of 
N2H+ and the J=l-0 and 3-2 lines of N2D+. Figure [12] 
shows a 5 points map of the J=l-0 transition of N2D+ 
(the reference position is given in table [Tj). Figure [2] com- 
pares the J=l-0 radial integrated intensity of the mo del 
to observational data reported by ICaselli et al.l l|2002f ) . 

In order to simultaneously reproduce the J=l-0 and 
3-2 lines of N2D"'", a step in the isotopologues abun- 
dance ratio has been introduced, this ratio varying from 



0.6lo 3 in the inner 10" to ' 



0.25_o;o2 at greater radii. 
Note that the extent of the enriched region favors the 
intensity in the J=3-2 line compared to the J=l-0 line 

(see discussion of L183). 

L63 has been mapped in NH3 ijBenson &: Mverslll989[ ) , 
N2H+ and CO (Sncll 19 8ll) as well as at 800 Mm and 1.3 
mm continuum emission ( Ward-Thompson et al.l [19941 . 
I1999I I. The intensity peak of both the continuum emis- 
sion and N2H+ maps are found at the same position. 
The central density and break radius obtained in this 
work (no = 7.2 10^ cm~^, ro = 0.018 pc) are in excel- 
lent agreement with the corres ponding values derived by 
IWard-Thompson et al.l (|l999h from the analysis of 1.3 
mm emission (no = 7.9 10^ cm~'^, ro = 0.017 pc). In 
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compa rison to the results derived by iBenson fc MversI 
l|l989t ) from NH3 observations, the profile determined in 
this work corresponds to a more massive cloud (within 
2.4', 8 M© from NH3 compared to 16.5 M© from N2H+). 
Such a discrepancy seems to be present all sources we 
have modeled. 

For this source, we have introduced a radial variation 
of the temperature, going from 8 K in the inner part to 15 
K in the outer envelope. This is consistent with previous 
works where the temperature h a.s been estimated to be 
8 K in the inner 40" f continuu m. IWard-Thomp son et all 
I|l994l l). 9.7 K inside 2.4' fNH^ iBenson fc Mve rs (1989)), 
and 15 K within 8' (CO. ISnelll l|l98lh ). Such an increase 
of the temperature in the outer regions of clouds is ex- 
pected due to heating of the gas and dust by the inter- 
stellar radiation field. 

3.1.7. L43 

Figure m shows the J=l-0, 2-1 and 3-2 Hues of N2H+ 
and the J=3-2 line of N2D"'", and Figure [2] presents 
the J=l-0 radial integrated intensity from the model 
compa red to the observations reported by Caselli et al.1 

Hooa). 

The dense core of L43 is composed of one main conden- 
sation with a filamentary extent at the South-East. This 
core is associated with the T-ta uri star RNQ91 located 
1.5' W est. As pointed out by Ward-Thompson et al.1 
l|l999t ). the initial cloud might be forming multiple stars, 
RN091 being the first to have appeared. From 1.3 
mm continuum emission, IWard-Thompson et al] l|l999l l 
found an average central volume density of ~ 2 10^ cm~ 
within a radius of ro=0.018 pc, at which a steeper slope 
is reached. The density profile we have determined is less 
centrally peaked (i.e ro = 0.022 pc, no = 9.2 10^ cm^^) 
but leads to the same estimate of the enclosed mass in 
the cloud: we obtain masses of 1.1 and 8.2 Mq within, 
respectively, 22" and 51". This latter radius corresponds 
to the geometri c average of the FWHM raajor and minor 
axis quoted by IWard-Thompson et al.l ()l999( ) at a dis- 
tance of 160 pc. Our results are in excellent agreement 
with the estimations of 1.2 and 8.2 Mq of the latter work. 

3.1.8. L1489 

Figure [H] presents the J=l-0, 2-1 and 3-2 Hues of 
N2H+ and the J=3-2 line of N2D+. Figure H shows the 
J=l-0 radial integrate d intensity conipared to observa- 
tional data obtained by ICaseUi et al.l l|2002l l . Our obser- 
vation of the J=l-0 line shows asymmetrical line wings 
with an enhanced red wing that we failed to reproduce 
in our models. Thus, the integrated intensity shown on 
Figure [2] is underestimated in the model by comparison 
to observational data. 

This source belongs to the Taurus complex and is as- 
sociated with the low mass YSO, IRAS 04016+2610 sit- 
uated ~ 1.2' to the West. The NH3 map shows one main 
condensation of radius ~ 0.07 pc with a peak intensity 
shifted ~ 1' North of the N2H+ peak. iBenson fc MversI 
l|l989l l estimated the temperature to be 9.5-10 K and the 
main condensation mass to be around 1.6-2.1 Mq, which 
is lower than our estimate of 6.5 Mq. On th e other hand, 
observ ations of HCO"'' Hues performed by lOnishi et al.1 
l|200^ have shown that two condensations are lying in 
the NH3 emitting region. These clumps have a low spa- 



tial extent (r = 0.018 and 0.016 pc ) and correspond to 
dense cores with masses ~ 1 Mq. 

3.1.9. L1251C 

Figure [H shows the J=l-0, 2-1 and 3-2 Hues of N2H+ 
and the J=3-2 line of N2D+. 

L1251C is a dense core belonging to a region of active 
stellar formation in the Cepheus complex. Large scale 
studies of t his reg ion are reported for ^"^CO and C^^O by 
ISato et al.1 l|l992h and for NH3 bv lToth et aH (|l996l l. As 
for the other clouds, we have determined a higher density 
from our N2H"'' data compared to that obtained from 
analysis based on other molecular species. This results in 
a larger mass estimate for t he cloud. In prev ious studies, 
masses of 5 Mq within 62" llToth et al.lll996l l and 56 Mq 
within 300" l|Sato et al.lll992l ) have been obtained. In this 
work, we derive masses of 13 and 87 Mq respectively. 

3.2. Discussion 

For all sources where the J=l-0 Hue of N2D+ was ob- 
served, we have been able to obtain models which re- 
produce the observed hyperfine intensity ratios. On the 
other hand, it is often more difficult to obtain such a 
good agreement when trying to model the J=l-0 hy- 
perfine line intensities of N2H"'". For this molecule, a 
number of studies have reported intensity anomalies for 
the J =l-0 hyperfine transitions (see e.g. ICaseUi et al.1 
Il995( ). Between the two isotopologues, the main differ- 
ence is due to optical thickness. Indeed, when comparing 
the observations of the J=2-l fines of N2H"'" and N2D"'", 
we see that, for most sources, we observe self absorption 
features in the J=2-l main central component of N2H"'" 
while not for N2D"'". Thus, this fact supports the idea 
that hyperfine intensity anomalies are due to radiative 
effects such as scattering in the low density regions of 
the clouds. Hence, the lower opacity of the J=l-0 line 
of N2D+ compared to N2H+ explains why such intensity 
anomalies are not observed for the deuterated isotopo- 
logue. 

For most sources, we find a better agreement between 
models and observations using a N2H+ abundance which 
varies with radius. The N2H"'' abundance is enhanced 
in the central region in all models, except for L1517B 
and for molecular peak C in L183. For L183, the best 
agreement with observations is found by introducing an 
increase of the N2H"'" abundance in the outer envelope. 
This could indicate a real trend for this source but, on 
the other hand, it must be confirmed with a more precise 
description of the morphology of the cloud, as previously 
discussed. Note that these two clouds are the less mas- 
sive in our sample (see table [l]). For all the other clouds 
for which we have obtained enclosed masses greater than 
5 Mq within 2', the N2H"'' abundance decreases in the 
outermost regions. These different types of abundance 
profiles were already noticed for some clouds in our sam- 
ple (see sect. [3]) and this refiects the specificity of the 
chemistry in each source. 

Compared to previous studies, we have derived lower 
N2H+ abundances despite the fact that the coUisional 
rates used in the present work are lower than the rates 
of HCO"*" - H2 currently used to interpret N2H"'" spec- 
tra. Models with high density and low abundance are 
the best to reproduce the J=2-l line which suffers from 
self-absorption. Moreover, in the sample of clouds we 



8 



studied, we find a weak departure in the abundances (less 
than a factor of 2) to the average N2H+ abundance of 1.5 
10~^°. Note that for the same sample of cl ouds, the av- 
erage value for the abundance reported in ICaselli et al.l 
(I2OO2!) is 4 10"^". In all these clouds, a decrease of the 
temperature in the inner region below 10 K enables to re- 
produce best the observations since it increases the opac- 
ity of the J=l-0 line. On the other hand, an increase 
of the temperature with radius allows to minimize self 
absorption effects in the J=2-l line. 

In our sample, several clouds have already been ana- 
lyzed using continuum observations (L63, L43, TMC-2, 
L183, L1498 and L1489). These studies lead to esti- 
mates for the mean volume density which are in good 
agreement with those derived in this work. Nevertheless, 
when using other molecules, the differences in density 
estimates are more important. Part of the discrepancies 
may arise from the different spatial distribution within 
the clouds of each molecule. Some of them, such as CO 
or CS, are known to undergo depletion in the highest 
density regions and an analysis which assumes a uniform 
density would lead to a lower value for the cloud den- 
sity. For the studies based on NH3, the differences are 
yet unclear since both molecules seem to be present in 
the same regions of the clouds, as shown by the corre- 
sponding emission maps. For the compar ison with the 
estimates done in iBenson fc Myer j (|l989h , a source of 
discrepancy could be the larger HPBW (~ 80") of the 
Haystack antenna used in the latter work. Another possi- 
bility explaining the lower density estimates derived from 
NH3 lines may come from the low critical densities of 
the (J,K) = (1,1) and (2,2) inversion lines, i.e. Uc ^ 10^ 
cm~^. Since for the densities of the inner regions of the 
dense cores, these fines are thermafized, an accurate de- 
termination of the kinetic temperature and NH3 column 
density is stfil possible while the fines are rather insen- 
sitive to the density. Moreover, the densities estimated 
from N2H"'" and millimeter continuum observations lead 
to NH3 l ine intensities cons istent with the observations 
(see e.g. iTafalla et al.l[20Q4f ). This points out the low 
sensitivity of the NH3 inversion lines to dense gas, and 
the difficulty to assess the dense core gas density uniquely 
from NH3 observations. Finafiy, we obtain that the mass 
derived for the clouds are in general in good agreement 
with the estimates obtained from mm or sub-mm emis- 
sion analysis and higher than the mass derived from the 
observations of other molecular species. 

Note, that the rate coefficients used in this work 
are calcu l ated f or collisions with He. As discussed in 
iMonteirol l|l985l l for the case of HCO"*" (isoelectronic of 
N2H+), the larger polarizability of H2 compared to He 
induces strong variations in the corresponding potential 
energy surfaces. It entails that the cross sections ob- 
tained for HCO+-H2 are greater by a factor 2-3 com- 
pared to the HCO"'"-He ones. In the present work, the 
collisional rate coefficients are corrected for the difference 
in the reduced masses of H2 and He in the Boltzmann 
average. This is a crude approximation which assumes 
that cross sections for collisions with He and H2 are sim- 
ilar. Thus, we expect to obtain significant differences, by 
analogy with the HCO"*" case, for the cofiisional rates of 
N2H"'" with H2. Therefore, the results obtained have to 
be regarded with caution and the error on the density 
and abundance estimates could be as large as a factor 2. 



The models presented for L1517B and L1498 can serve 
as a benchmark for the typical variations that we expect 
using different collisional r ate coefficients. Com pared to 
the abundances reported in lTafalla et al.l l|2004f ). the val- 
ues we have derived in this work are within 50-70%. 

From N2D+ fines, we have derived values for the D/H 
enrichment in each cloud which are in reasonable agree- 
ment with previous determinations, despite the fact that 
the methods used are often different. In addition, for 
three of the studied clouds (L63, TMC2 and the two cores 
in L183), we find that the D/H ratio takes high values 
(0.5-0.7) in the inner part of the cloud and decreases in 
the lower density outer regions. Such characteristics are 
expected from theoretical models which c ouple the dy- 
namics and chemistry of molecular clouds (jAikawa et al.l 
[2005; Robe rts et al, 2QQ3i ). Moreover, these models pre- 
dict that a quantitative determination of the deuterium 
enrichment could probe the evolutionary stage of the col- 
lapse. 

4. CHEMICAL ANALYSIS 

The present analysis of the line profiles of N2H"'" and 
N2D"'" has allowed to derive the profile and temperature 
dependences together with the overall deuterium frac- 
tionation ratio in the various selected pre-stefiar cores. 
Looking at tabled one may notice that the deduced to- 
tal column density of N2H+ varies within a factor of 2 
whereas the deuterated counterpart abundance may vary 
by factors larger than 6, which leads to fractionation val- 
ues between 0.07 and 0.5. 

A detailed chemical discussion on the deuterium frac- 
tionation derived from the present observations is be- 
yond the scope of the present paper and we only want 
to derive some general trends from the observations. 
The presence of molecular ions in interstellar envi- 
ronments is readily explained by a succession of ion 
molecule reactions initiated by cosmic ray (CR) ion- 
ization of molecular hydrogen : H2 + CR — > Hj; H2 
+ B.2 Hj^ + H. The H^ molecular ion is an ef- 
ficient proton donor and reacts with saturated stable 
molecules such as CO, N2, HD to produce HCO+, N2H"'", 
H2D''". Deuterium enhancement follows when the tem- 
perature is low enough so that + HD deuterium ex- 
change reaction may proc eed on the exo thermic path- 
way as first proposed by IWatsorJ l|l974l l. The H2D+ 
thus produced can further transfer its deuteron to the 
same neutral molecules, producing DCO"*", N2D"'", at 
a rate 3 times smaller than in the reaction involving 
B.^ , as 2/3 of the reaction rate coefficient will lead 
also to HCO"^ and N2H"'", if statistical arguments may 
apply. In the reaction of H2D+ with HD, D2H+ is 
formed, and a further reaction with HD then leads to 
the completely deuterated ion D;^. These p rocesses have 
received attention both in astro phvsics ([Roueff et ap 
20001 [Rob erts fc Mifia^ l2000l : [Roberts et all l20Q3l : 
Flower et al . 2004) i an d in chemical physic s stud - 
ies (iGerfich et all l2002t iGerlich fc Schlemmerl l2002h . 
iRamanlal fc TennvsonI ( 20041 ) have summarized the 
exothermicities involved in the possible reactions involv- 
ing deuterated substitutes of H3 and H2 , which are typ- 
icafiy between 150K and 220K. Temperatures as those 
derived in Section[3]are thus low enough to inhibit the oc- 
currence of the corresponding reverse reactions and suit- 
able for deuterium enhancement to occur. The confirmed 
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detections of H2D" '" in the pre- s tellar core LDN 1544 
ijCaselli et alJ l2003l : IVastel et alJ l2006l l. as well as the 
recen t detection of D2H+ in LDN 1689N ijVastel et all 
l2004f ) , have provided strong observational support for 
this theory. It is important to notice that the isotopo- 
logue ions thus formed (N2H+, N2D+ on the one hand, 
HCO"'", DCO"'" on the other hand) are then destroyed by 
the same chemical reactions, i.e., dissociative recombina- 
tion, reactions with the other abundant neutral species, 
possible charge transfer reactions with metals, recombi- 
nation on grains, ... In the absence of specific measure- 
ments, one assumes that the rate coefficients involving 
hydrogenated and deuterated molecular ions are equal. 
The steady state fractionation ratios N2D''" / N2H''", 
DCO"'' / HCO"*" are then directly equal to the ratio of 
the formation probabilities of N2H"'" and N2D~'' on the 
one hand and HCO"*" and DCO"*" on the other hand. 

In the case of N2H"'", the main formation route is usu- 
ally N2 + Hj^, whereas for N2D+, the formation channels 
involve N2 + H2D+, N2 + D2H+ and N2 + D;^ . N2H+ 
may also be produced via N2 + H2D+ and N2 + 0211+ if 
these ions become abundant. At steady state, the abun- 
dance ratio is given by : 

X(N2D+) ^ X(H2D+) + 2 X(HD+) + 3 X(D+) 
X(N2H+) ~ X(HDj) + 2X(H2D+) + 3X(H;[) ' 

where X(x) represents the abundance or the fractional 
abundance of a particular species x. This relation holds 
if the overall reaction rate coefficients of N2 with the 
various deuterated isotopologues of H;^ are equal and if 
statistical equilibrium determines the branching ratios. 

High fractionation ratios such as those found in 
the studied environments are then directly correlated 
to high deuteration fractions of the deuterated 
ions which can occur in highly depleted environments 
llRoberts et~an! 2004: Walmslev eraLll20Q4l : iFlower et all 
l2004t [Roueff et al.ii2Q05h . In addition to these reactions, 
N2D"'" (DCO"*") may also be formed in the deuteron ex- 
change reaction between D atoms and N2H"'" (HCO"*"). 

Density and temperatures are derived from the anal- 
ysis of the line profiles as shown in Table [2j The de- 
duced fractional abundances of N2H'^, which are of the 
order of 10~^° compared to H2, imply a first constraint 
on the depletions, i.e. the available abundances of the 
elements in the gas phase. The deuterium fractionation 
ratio brings an additional limit on these values and also 
on other physical parameters such as the cosmic ray ion- 
ization rate {Q which is directly related to the ionization 
fraction. 

Steady state model calculations allow to span rapidly 
the parameters space and to obtain the main phys- 
ical and chemical characteristics of the environment. 
We have use d an u pdated chemical model described in 
iRoueff et al.1 l|2005l l. where we include multiple deuter- 
ated species containing up to 5 deuterium atoms (CD^) 
and the most recent chemical reaction rate coefficients. 
The gas phase chemical network includes 210 species con- 
taining H, D, He, C, N, 0, S and a typical metal un- 
dergoing charge transfer reactions which are connected 
through about 3000 chemical reactions. Particular care 
is paid to dissociative recombination branching ratios^ 

^ In the current chemical analysis, the branching ratios of dis- 



which modify considerably the chemical composition of 
the gas. We mimic the various mechanisms involved 
in gas-grain interactions (accretion, thermal desorption, 
cosmic-ray induced desorption, etc ...) by introducing a 
dependence of the available gas phase elemental abun- 
dances with density. The variation laws assumed for the 
gas phase elemental abundances with density are arbi- 
trary and account for the decrease of gas phase atomic 
carbon, oxygen, nitrogen and sulfur with increasing den- 
sity. The assumed laws are: 

nc/uH = 1.40 10"^ [1 - exp (-4 10* /uh)] 

riN/riH = 7.94 10"^ [l - exp (-3 10* /hh)] 

no/nn ^3.5010-* [l - exp (-4 lOV^ff)] 

ns/uH = 1.85 10"^ [1 - exp (-4 10* /n„)] 

where rix stands for the elemental volume density of 
species x. At low densities, the values correspond to the 
abundances observed in diffuse and translucent environ- 
ments. The carbon to oxygen ratio is fixed throughout 
the density variations. For H2 densities in the range be- 
tween 10^ and 10^ cm~^, the elemental abundances of 
C, N, O and S are in the range of several 10~^ - 10~^, 
corresponding to depletions of 10 - 100. 

We solve the charge balance on the grains and allow 
the atomic ions to neutrafize on the mostly negatively 
charged grains. These processes are important for the 
overall ionization fraction of the gas and lead to very low 
ionization fraction of the order of 10~^ for densities larger 
than 10^ cm~'^. S uch low value s are c ompatible with 
the derived value bv lCaselli et al.l (|2002br ) from observed 
deuterium fractionation in dense cores. We show in Fig. 
[161 both the fractional abundance of N2H+ in units of 
10~^° and the N2D"'" / N2H"'" ratio for densities ranging 
from 10^ to 10^ cm~ and a temperature of lOK. Two 
( values are used to probe the role of this important 
physical parameter. 

We see that the obtained results span very nicely the 
values derived from the observations. To model the 
N2H+ emission, it has been assumed that the density 
structure was consisting on a power law. We thus derived 
that observed spectra were consistent with an increase 
of the temperature with radius, and that, depending on 
the source, the N2H+ abundance was whether constant 
or decreasing outward. Considering the two ( values as- 
sumed we find that the trend derived from the modeling 
is qualitatively reproduced with C = 1 10"^'' s~^. On the 
other hand, for a higher ( value, the abundance would 
tend to increase with radius. 

The models done in section [3] show that the 
N2D+/N2H+ ratio could reach values as high as 0.5 in 

sociat ive recombination of N2H+ are taken from IGeppert et alj 
II2004I V During the revision of the manuscript, it has been brought 
to our attention that in in a recent experiment, [Molek et al. (2003) 
have reanalysed these branching ratios and o btained a major chan- 
nel of reaction towards N2, opposite to what IGe ppert et al. ( 2003) 
had found with the storage ring experiment. For some test cases, 
we have introduced these new branching ratios, by keeping the 
same value of the total destruction rate of N2H+ and N2D+ and 
we have found that the fractional abundances of N2H+, N2D+ as 
well as the fractionation ratio remain very similar to the one dis- 
played on figure 1161 However, substantial changes are found for 
hydrogenated nitrogen molecules and their isotopologues such as 
NH, NH2, ... 
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the innermost regions of the clouds. FigurefTBl shows that 
chemistry predict a N2D+/N2H+ ratio greater than 0.1 
for densities higher than 2-3 10^ cm~'^, with maximum 
values of 0.2-0.3. While the maximum values derived 
from the chemical modeling and from observations are 
within a factor 2, the trend derived for the ratio, in both 
cases, are in good agreement. Actually, for four sources 
(L63, TMC2, and the two dense cores of L183), we ob- 
tained from observations that the ratio should decrease 
rapidly with density, a trend reproduced in the chemical 
modeling. 

Note that both absolute abundances and abundances 
ratio have been found to be sensitive to the adopted 
depletion law. In contrast. Figure [16] shows that the 
N2D"'"/N2H'^ ratio is rather insensitive to variations of 
the temperature or cosmic ray constant, while the N2H+ 
and N2D"'" absolute abundances are strongly changed. 
We consider here only gas phase processes so that the 
role of induced desorption by cosmic rays is assumed to 
be constant within the variation by a factor of 5 of the 
cosmic ray ionization rate value. A variation of C changes 
both the abundances of HJ^ isotopologues and of the elec- 
trons, which are involved respectively in the main N2H+ 
and N2D+ formation and destruction routes. The larger 
efficiency of formation from HJ^ isotopologues entails that 
the absolute abundances increase with (. Moreover, the 
similarity of the formation and destruction channels of 
N2H"'" and N2D"'" lead to nearly equal variations, both 
with the temperature and (, which explains why the iso- 
topologues abundances ratio remain identical. Finally, 
we point out that the N2D"'"/N2H"'" ratio could serve as 
a probe of the way the C,N,0 and S containing species 
deplete onto grains, since other parameters such as tem- 
perature and the cosmic ray ionization rate just influence 
the isotopologues absolute abundances similarly. 

Figure flGl shows the N2D+/N2H+ ratio derived from 
the full chemical network in comparison to the ratio ex- 
pected from statistical considerations (i.e. equation [1]) . 
We find a reasonable agreement between the two ratio, 
the largest differences (for T = 15K and C = 5 10"^'' 
s~^ ) being of the order of 60%. The largest values in 
the full modeling are due to an increased N2D"'" forma- 
tion through the deuterium exchange reaction between 
D atoms and N2H+. For some specific conditions, this 
formation route is found efficient. 

5. CONCLUSIONS 

We have derived estimates of the temperature, den- 
sity, N2H"'" and N2D"'" abundances in a sample of cold 
dark clouds, by using a radiative transfer modeling which 
enables to interpret the hyperfine transitions of these 
molecules. The main conclusions are : 

1. Compared to previous studies considering N2H+, 
we generally derive higher densities due to the in- 
clusion in the models of the J=2-l line. For the 
conditions prevailing in dark clouds, i.e. Tk ~ 10 
K, this line is the optically thickest and, in order 
to prevent strong self absorption effects which are 
generally not observed, models with high density 
and low abundance are preferred. Moreover, an 
outward increase of the temperature enables to re- 
duce the self absorption in this line. 

2. In order to reproduce the J=l-0 hyperfine transi- 



tions, the total opacity of this line has to be larger 
than r ~ 10. A way to increase the opacity and 
to still reproduce the other rotational lines is to 
decrease the temperature in the inner core below 
T = 10 K. For most sources, we encounter a good 
agreement with a temperature around T = 8 K in 
the central region of the clouds. 

3. We analyzed the density structure of the clouds 
and we took into account previous studies where 
the J=l-0 integrated intensity, as a function of the 
position on the cloud, were reported. The central 
average densities and the radii of the inner fiat re- 
gions, derived using N2H"'", are in good agreement 
with the equivalent parameters obtained from sub- 
mm and mm continuum observations. Compared 
to other studies based on other molecular species, it 
seems that N2H"'" is the only one which allows such 
a good agreement, since other molecules provide 
systematically lower H2 volume densities. For most 
molecules, it still might be possible to conciliate 
these estimates by introducing the radial depen- 
dence in the molecular abundances which are pre- 
dicted by theoretical studies of dark clouds chem- 
istry. 

4. We derive X(N2D+)/X(N2H+) which are in qual- 
itatively good agreement with previous studies. 
Moreover, for two of the studied clouds (TMC2 and 
L183 ) where we have observed various positions in 
the cloud, we find that the observed spectra are 
best reproduced when increasing the D/H ratio in 
the inner dense regions. We note that for these 2 
objects a single position analysis lead to the deriva- 
tion of a constant ratio throughout the cloud, which 
was smaller than the central value obtained in the 
multi position analysis. Thus, for the other objects, 
the ratio we have derived has to be considered as 
an average value and a central enhancement of the 
D/H ratio cannot be ruled out. 

Considering studies that deal with the chemi- 
cal evolution of protostellar clouds, our results 
are in qualitative and quantitative agreement 
with the expected t rend of the two isotopologues 
(jAikawa et al.l [20051 ). Recently, the importance of 
the multideuterated species of Hj" in the chemica l 
network was pointed out by [Roberts et al.l ([2003r ) 
in order to explain the high deuterium fraction- 
ation observed in molecular clouds. This work 
shows that for various molecular species we should 
observe D/H ratios approaching unity. Moreover, 
they show that the deuterium fractionation is par- 
ticularly efficient in the case of N2H"'". A detailed 
time dependent study consid ering both dynamics 
and chemistry was reported in lAikawa et al.1 (|2005D 
for the stages prior to star formation. It confirms 
that deuterium fractionation can reach high values 
for N2H+ and that the abundance ratio takes larger 
values in the innermost regions of the clouds. From 
this latter study, it appears that the fractionation 
is expected to increase with time when the central 
density of the inner region increases. Thus, the 
determination of the deuterium enrichment could 
serve as a tool to probe the stage of evolution prior 



11 



to the formation of protostars. 
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Appendix: Radiative Transfer Models for N2D"'" 

The rotational and hyperfine c oupHn g constants for 
N2D"'" are taken from iDore et al.l (|2004l ). We have used 
the dipole moment of N2H"'", i.e. fi = 3.4 D, to derive 
the line strengths of the N2 D"'" hyperfine Hues, according 
to equation (2) of iPaperl . For the hyperfine collisional 
rate coefficients, we have used the de-excita tion rates of 
N2H+ colliding with He ([Daniel et al.ll2005f ) and applied 
the detailed balance relationships to derive the excitation 
rate coefficients (see below). 

Figure [17] shows the excitation temperatures, opacities 
and brightness temperatures of the J=l-0, 2-1 and 3-2 
transitions of N2H+ and N2D+ obtained under the LVG 
approximation. The column density is the same for both 
species and equal to 10^^ cm~^ / (km s~^ pc~^). In order 
to derive the rotational Hue opacities we hav e to take into 
account the hyperfine structure (see PapeTll). Excitation 
tempera tures are obtained according to equation (7) of 
iPaper J . 

The difference in the excitation processes for the two 
isotopologues is a consequence of their specific rota- 
tional energy structure, and hence, different frequencies 
and Einstein coefficients for the rotational lines. For 
N2D"'", the frequencies of the rotational transitions are 
lower than those of the main isotopologue by a factor 
B(N2H+)/B(N2D+) = 1.21. Consequently, the Einstein 
coefficients are lower for N2D+ by a factor 1.8, which re- 
sults in lower critical densities if we assume that the de- 
excitation rate coefficients are identical for both species. 
This latter assumption has been checked by comparing 
the rotational rate coefficients of N2H+ and N2D''" col- 
liding with He (see table [5]). The N2D"'" rate coefficients 
have been computed by the MQLS CAT^ code u s ing th e 
potential energy surface given by iDaniel et al.1 l[2004l ). 
The reduced mass of the colliding system has been mod- 
ified with respect to the N2H'^ calculations. However, we 
have assumed that the mass center for both isotopologues 
were the same (this enables a straightforward determina- 
tion of rate coefficients for the deuterated species) . The 
results are given in Table [H We can see that for temper- 
atures ranging from 10 to 30 K, the differences between 
the rotational de-excitation rate coefficients of the two 
isotopologues are < 20%. However, upwards rate coeffi- 
cients are obtained from detailed balance and they will 
be different for both isotopologues, being larger for the 

^ J. M. Hutson and S. Green, MOLSCAT computer code, version 
14 (1994), distributed by Collaborative Computational Project No. 
6 of the Engineering and Physical Sciences Research Council (UK). 



deuterated species. The difference is small, of the order 
of 10-20%, for AJ=+1 transitions, but becomes as large 
as a factor 4 for higher AJ. The lower Einstein A co- 
efficients and higher upwards rate coefficients for N2D+ 
with respect to N2H+ allow the high-J rotational lev- 
els to be more efficiently populated for the deuterated 
species. Hence, the effect introduced by the different 
rotational constants has to be taken into account when 
computing the level population for both isotopologues. 
Figure [TTI shows that for identical column densities (10^^ 
cm~^) the J=l-0 and 3-2 lines of N2D"'" have, respec- 
tively, lower and higher opacities than those of N2H+. 
This behavior applies to the whole range of densities ex- 
plored in our calculations. Nevertheless, the change in 
the critical density is only apparent for the J=l-0 line 
(see Tex panel of Figure fT7[) . The effect is less evident for 
the other transitions because rotational excitation tem- 
peratures are obtai ned frorn the average of the individual 
hyperfine Tex (cl. iPaper J ). 

The expected variation of critical densities is more ob- 
vious when considering HCO"*", H^^CO"'' and DCO+. 
Figure [17] (right panels) shows that the opacities for 
each transition vary similarly than N2H+. More- 
over, on the Tex panel of Figure [17] we see that the 
critical densities are lower for DCO"*" and H^'^CO"'" 
than for the main isotopologue. It is worth noting 
that the intensity ratios Tb(DCO+)/Tb(HCO+) and 
Tb(N2D+)/Ts(N2H+), depend on the volume density 
for 10'' < n(H2) < 10° cm~^ (see bottom panels of Fig- 
ure [TT]). This behavior is found even for higher column 
densities. If the column densities for the deuterated and 
main species are not the same, then the plots in the bot- 
tom panels of Figure [17] should be scaled by the assumed 
column density ratio. 

The effect discussed above for singly deuterated lin- 
ear molecules could be even more important for multiply 
deuterated species. In order to evaluate qualitatively this 
effect we have considered the excitation of two asymmet- 
ric top molecules, 0-H2CO and P-D2CO (note that the 
quantum numbers for ortho-formaldeyde and its dou- 
bly para-deuterated isotopoloque are the same). Fig- 
ure [18] shows the excitation temperature in the optically 
thin case (column density of 10^^ cm~^ for both isotopo- 
logues) of different rotational transitions of these species 
as well as the corresponding line intensity ratios. The 
difference in the Einstein coefficients for P-D2CO and o- 
H2CO, and the larger upwards excitation rates for the 
former, produce an important variation in the excitation 
conditions of these two molecules. The line intensity ra- 
tio varies from 1.7 for the transition 2i2-lii to 3.4 for 
the 4i4-3i3 one. We stress that this difference in the ex- 
citation conditions will translate into different emitting 
volumes for each species. Hence, the determination of 
the deuteration degree can not be directly related to the 
column density ratios of the main and deuterated iso- 
topologues. The cloud density and temperature struc- 
ture has to be used, together with non-local radiative 
transfer codes, in order to derive from observations the 
deuteration degree through the cloud. 

The effect will be even more important for triply 
deuterated species like CD3OH or CD3CN for which the 
rotational constant A will change b y a factor :;^2 .0 and 
B and C by a factor ~1.5 (see e.g. IWalsh et al.lll998.) . 
Thus, the energies of the CD3OH levels are lower by a 
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factor ~ 2 compared to CH3OH. Hence, from detailed 
balance relationships, upward collisional rates will be 
found very different for both isotopologues. The global 
effect in the determination of the deuteration enrichment 
will also depend on the source physical structure. While 
the high observed beam averaged abundance ratios for 
CH3OH and its doubly and triply deuterated isotopo- 
logues can not be ruled out by these excitation effects (see 



e.g. iParise et al.ll2002l [2003 l. an exact determination of 
these ratios should take into account the different cloud 
emitting volumes for these species. The implication for 
chemical models is obvious since a direct comparison of 
the derived column densities will always produce higher 
deuteration factors than the real ones. Thus, detailed ra- 
diative transfer models have to be carried out to derive 
correct values for the deuteration enrichment. 
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Source 


«(J2000) 


(5(J2000) 


I(K km s"^) 


rf Dc) 


D(pc) 


L1489 


04 04 49.0 


26 18 42 


1.8 


0.039 


140 


L1498 


04 10 51.4 


25 09 58 


1.5 


0.051 


140 


TMC2 


04 32 46.8 


24 25 35 


2.5 


0.061 


140 


TMC1(NH3)^ 


04 41 21.3 


25 48 07 


2.3 


0.040 


140 


L1517B 


04 55 18.8 


30 38 04 


1.4 


0.050 


140 


L183(C) 


15 54 08.7 


-02 52 07 


2.5 


0.045 


110 


L183(N) 


15 54 09.2 


-02 49 39 


1.8 


0.045 


110 


L43 


16 34 35.0 


-15 46 36 


3.4 


0.059 


160 


L63 


16 50 15.5 


-18 06 26 


2.0 


0.057 


160 


L1251C 


22 35 53.6 


75 18 55 


1.7 


0.069 


200 



TABLE 1 

Right ascension and declination of the observed sources, assumed distance (D), integrated intensity of the J— 1-0 line 

(I) AND half power RADIUS OF THE MODELS (r), AS OBTAINED BY CONVOLUTION WITH THE FCRAO ANTENNA BEAM, FOR 
COMPARISON WITH THE OBSERVED VALUES GIVEN IN ICaSELLI ET AlI ()2002f l. 

^ For this source, the half power radius and the central intensity are calculated considering only the main emission 

component. 



Source 


M"] 


no/10= [cm-^] 




Ti[K] 


Xi/10-1" 


V2["] 


T2[K] 


X2/IO-"' 


■■3["] 


T3[K] 


X3/10-" 


L1251C 


29+1 


°-°-0.5 


21 


7 + 0-2 
' -0.4 


i.5is-.;3 


150 


,(,+1.1 
J- '-'-2. 3 


o.2t°:; 








L43 


29+; 


q 2+°-* 


30 


„ + 0.3 
°-0.9 


0.812.; 


60 


,,+1.8 

J- .1-2.0 


o.4t°:; 


150 


1411 


3+"-^ 
'J'-'-O.l 


TMC2 


60+t 


o r+O-l 


40 


+ 0.2 
''-0.6 


1 1+° 1 


200 


jO+1.8 
'^'^-1.4 


2+°-i 








L63 


23+; 


'•^-0.5 


20 


+ 0.3 
°-0.6 


1-41;!.;. 


45 


+ 0.8 
°-2.0 


o.4t°-.; 


150 


15 + 1 




L1489 


■24t\ 


_ O + 0-3 
5.8_(, 4 


35 


+ 0.6 
°-0.7 


0.8l°-i 


65 


q + 2.8 
^-1.6 


0.8l°:t 


100 


7 


1+"-^ 


L1498 


70+1 


04+0 "3 


70 


+ 0.3 
°-0.2 


2 5 + °-l 


200 




3 + °-i 
^■•'-0.1 








TMC1(NH3) 


30+; 


'■°-o.i 


35 


O + 0.4 
°-1.6 


, 0+0.1 

^■■'-0.3 


65 


,2 + 1-8 
J-^-1.6 


, + 0.2 
.'■■'-0.1 


175 




1+"'^ 


L1517B (1) 


30+; 


1 9+°-i 


18 






180 




1 6+°-i 
^•°-o.i 








L1517B (2) 


35+; 




18 




2 5+"-'' 


180 




e;+0.2 

^■■'-o.i 








L183(C) 


19+; 


R + O-2 


20 


O + 0.4 
''-0.6 


1 4+°-i 


40 


q + 0.6 
^-1.6 


1 4+°-l 
.'-■*-0.6 


225 


"-3 


2 1+"-^ 


L183(N) 




1 4+°-i 


80 


q+i.i 


2 Q+0.6 


225 


9l^ 


2 0+*" 









TABLE 2 

Parameters obtained from the modeling of our observations of N2H+ lines for the 10 clouds studied. The density profile is described as n(H2) = no for r < ro 

AND n(H2) = No(ro/r)^ for R > RQ. DEPENDING ON THE SOURCE, WE HAVE INTRODUCED 2 OR 3 COMPONENTS TO DESCRIBE THE CLOUD: THE RADIUS OF THE I*** COMPONENT IS 
GIVEN AS Rj AND FOR RADII T IN THE RANGE rj — l < T < Ti, THE N2H+ ABUNDANCE AND GAS TEMPERATURE ARE RESPECTIVELY GIVEN AS X, AND Tj (tHE FIRST ZONE WITH Tl AND 
Xl, CORRESPONDS TO THE RADIUS RANGE < r < ri"). FOR TWO CLOUDS, I.E. L1489 AND TMC1(NH3), THE TEMPERATURE OF THE 3'''* ZONE IS POORLY CONSTRAINED BY THE 
AVAILABLE OBSERVATIONS. ThE LARGE ERROR BARS FOUND MAKE THE VALUE DERIVED NON SIGNIFICANT AND THIS PARAMETER IS THUS NOT INDICATED. 



Source 


M(r<2') [Mq] 


rn 


N(N2H+)/10i^ [cm-=] 


N(N2D + )/10i^ [cm-=] 


R 


L12B1C 




n M+o-''^ 
J^*-o.ii 


12.6+1-1 


l.D_ J 3 





, Q + 0.06 
-'-■'-0.11 


L43 


26. 2t;;^ 


(1 lo+o-"* 


14.111.^3 


1 7+0.6 





, „+0.06 

-'^-o.oe 


TMC2 


22. 7+^7 


' ■'-0.38 




^■^-2.1 





rq+0.52 
^■'-0.28 


L63 


i3.4t;;i 




12.7l«:J 


c 1+2.7 
^■J--l,6 





^„+0.33 
^'^-0.14 


L1489 


- O + 0-6 

'■''-o.e 


"■ ' -0.11 


7.7t1i 


1 n + 0.5 

^^■•'-o.s 





1 7+0.12 
^ ' -0.11 


L1498 


7 6 + ^-° 

'■0-0.8 


"■" ' -0.02 


' ■'^-0.2 







07+0.02 

" ' -0.02 


TMC1(NH3) 


'■-'-0.4 


"■ ^^"^-0.03 


„ n+O.S 
"■"-i.o 







,„+0.09 
^^^-0.04 


L1517B (1) 




"■^'^-0.02 


6.3_o 4 


0.8l°:l 





, 0+0.04 
-^^-0.03 


L1517B (2) 




"■ ^^"^-0.02 










L183(C) 


■'■°-0.3 


"■^■'-0.12 




4 2+1'^ 
^■^-1.1 





,c, + 0.16 
''■'-0.11 


L183(N) 


- + 0-4 
' -0.4 


„ »„+1.00d 
U.7b_o 20 




1 4+1-* 

^^■^-0.2 





-1 c+o.ie 
-'^-o.oe 



TABLE 3 

For the clouds with model's parameters given in table[21 the mass enclosed within a radius of 2' is given (in Mq) as well as the abundance ratio no = 
X(N2D+)/X(N2H+). The column densities derived for the positions reported in table[T]are quoted and their ratio, R, is given in the last column. 

" THE quoted abundance RATIO IS FOR RADII BELOW 35" AND IS DECREASED TO 0.10 + q'5o FOR GREATER RADII. 
'' THE QUOTED ABUNDANCE RATIO IS FOR RADII BELOW 10" AND IS DECREASED TO 0.25lgj2 FOR GREATER RADII. 
THE QUOTED ABUNDANCE RATIO IS FOR RADII BELOW 40" AND IS DECREASED TO O.OS + gQg FOR GREATER RADII. 

THE QUOTED ABUNDANCE RATIO IS FOR RADII BELOW 10" AND IS DECREASED TO 0.05 + g q1 ^^R GREATER RADII. 
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N2H+ 




N2D+ 




Source 






LS R 


^ ^LSR 


Av 




(J=l-0) 


(.J=2-l) 


(J=3-2) 


(J=3-2) 


(km pc""'^) 


L1489 


6.74 


0.02 


0.07 


0.04 


-1.3 


TMC2 


6.17 


0.02 


0.01 


0.06 


-0.9 


T 18'i('C'l 


2.40 


0.01 


0.04 


0.04 


-0.2 


L183(N) 


2.41 




0.01 


0.06 


-0.0 


L43 


0.74 


0.03 


0.04 


0.01 


-0.4 


L63 


5.73 


0.02 


0.04 


0.03 


-0.0 


L1251C 


-4.74 


0.02 


0.02 


0.02 


-0.5 


L1498 


7.80 








-0.0 


TMC1-NH3 


5.98 








-0.8 


L1517B 


5.78 


0.04 






1.4 



TABLE 4 

Vz,SH OBTAINED FOR THE DIFFERENT TRANSITIONS OF N2H+ AND N2D+, BY MINIMIZING THE CHI SQUARE BETWEEN OBSERVED AND 
MODELED LINES. AlL THE VALUES ARE GIVEN IN COMPARISON TO THE VlSR OF THE J=l-0 LINE, I.E. VlSR = ViSi{(J=l-0) + 5 
^LSR- ^ V CORRESPONDS TO THE VELOCITY GRADIENT USED IN THE MODELS. NOTE THAT L1517B IS THE ONLY CLOUD WITH AN 

OUTWARD GRADIENT. 







T=10K 


T=20K 


T=30K 




j 


j' 




) -5 (%) 




-) 8 (%) 


R,^,/(N2D+) 5 (%) 


1 





103.39 


-10.6 


91.96 


-10.7 


88.51 


-9.4 


2 





58.20 


-7.9 


47.91 


-6.4 


42.09 


-4.0 


3 





35.04 


-4.7 


29.15 


-6.5 


24.98 


-9.3 


4 





19.81 


-16.5 


17.97 


-10.4 


16.44 


-9.4 


5 





19.33 


-13.0 


18.65 


-12.8 


18.01 


-12.2 


6 





21.78 


-2.7 


19.41 


-5.2 


17.84 


-4.9 


2 


1 


175.25 


0.0 


156.25 


-2.4 


145.50 


-3.3 


3 


1 


107.85 


5.7 


89.89 


-1.1 


78.69 


-3.7 


4 


1 


64.71 


-7.0 


59.22 


-7.6 


54.97 


-8.4 


5 


1 


61.24 


-13.9 


55.96 


-13.8 


52.23 


-13.0 


6 


1 


59.52 


1.6 


53.17 


-3.2 


48.58 


-4.6 


3 


2 


172.15 


-5.8 


155.22 


-5.7 


147.37 


-5.4 


4 


2 


125.74 


-11.7 


114.11 


-12.9 


106.30 


-12.8 


5 


2 


112.87 


-2.3 


100.39 


-4.8 


90.80 


-6.3 


6 


2 


79.06 


8.0 


73.04 


2.1 


67.56 


-0.9 


4 


3 


196.35 


-2.3 


182.35 


-3.4 


173.38 


-3.9 


5 


3 


127.69 


3.0 


119.03 


0.2 


111.16 


-2.1 


6 


3 


95.84 


4.2 


93.03 


2.8 


88.18 


0.5 


5 


4 


140.44 


-2.2 


141.22 


-0.5 


141.61 


-0.6 


6 


4 


110.55 


3.4 


108.73 


1.8 


104.69 


-0.3 


6 


5 


119.72 


-3.2 


122.05 


-0.7 


124.50 


0.0 



TABLE 5 

Rotational rate coefficients for N2D+ colliding with He, in units of 10"-*^^ cm* s"'^, and percentage difference with 

THE N2H + -HE de-excitation RATE COEFFICIENTS. 
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J=3-2 
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Fig. 1. — Observed 7=1-0, 2-1, and 3-2 transitions of N2H+ in the direction of selected sources. Ordinate is antenna temperature and 
abscissa is the Ylsr- For each transition the assigned frequency corresponds to the one of the strongest hyperfine component. 
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Fig. 2. — Observed rad ial integrated inten sity profiles of L1517B, L1498, L1489, L43 , L63 and TMC2 c ompared to the models. 
Observational data are from ICaselli et all ^2003 ) (obtained with the FCRAO antenna) and ITafalla et aP ||20Q4I ') for L1498 and L1517B 
(30m, Pico Veleta). For each source, the solid lines refer to the emerging intensities of our modeling (see text). For L1517B, the solid line 
correspond to the density profile used in this work and the dashed one to the profile used in ITafalla et alj 1)20041 '). 
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Fig. 3. — Observed and modeled J=l-0 and 2-1 transitions of N2H+ and N2D+ towards L1517B. T he solid line corresponds to the 
density profile used in this work and the dashed one corresponds to the profile used in lTafalla et alj 1)20041 '). 
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Fig. 4. — Observed and modeled J=l-0, 2-1 and 3-2 lines of N2H+ and J=2-l and 3-2 lines of N2D+ towards molecular peak C in 
L183. 
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Fig. 5. — Observed and modeled J=l-0 and 3-2 lines of N2H+ and J=l-0, 2-1 and 3-2 lines of N2D+ for the molecular peak N in 
L183. 
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Fig. 6. — Map and modeled lines of the J=2-l transition of N2D+ corresponding to molecular peak C in L183. The reference position 
of the map is reported in table [T] and corresponds to the N2H+ and N2D+ observations shown on figure ID On the right side of the main 
map, the map of integrated intensity is reported with contour levels going from 20% to 90% of the map peak intensity, by step of 10%. 
Note that the aspect ratio between right ascension and declination is respected on the right map and show the South-North extent of the 
cloud. 




^30 -20 -10 

right ascension 

Map and modeled lines of the J=l-0 transition of N2D+ corresponding to the molecular peak C in L183. 
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Fig. 8. — Observed and modeled lines for two positions toward TMC2. At the position of the N2H+ map peak are shown the J=l-0, 
2-1 and 3-2 lines of N2H+ and the J=3-2 line of N2D+. For the offset position, are shown the J=l-0 line of N2H+ and the 7=2-1 line 
of N2D+. 
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Fig. 9. — Models for the J=l-0, 2-1 and 3-2 line of N2H+ and for the J=3-2 line of N2D+ observed in TMC1(NH3). The upper right 
panel is a blow-up of the low frequency hyperfine triplet in the J—l-0 line, with the second emission component indicated. 
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Fig. 16. — N2H+ fractional abundances (solid lines, in unit of 10 ^'^) and ratio of abundances r£)=X(N2D+)/X(N2H+) (dashed lines), 



for a cosmic ray constant of C = 1 10' 



-17 



(upper panel) and C = 5 10 



-17 



-1 



(lower panel). For each panel, the results are shown for 



clouds at T = 7K (faint grey), T = lOK (grey) and T = 15K (black). The lower-right box of each panel compares the ratio ro to the ratio 
rlJ"* predicted using equation [T] Central densities (no) , N2H+ abundances, and abundance ratio (R) derived from models (see section [3j 
and listed in tables |2] and [3] are indicated in the lower panel, for comparison with the chemical models. 
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Fig. 17. — J=l-0, 2-1 and 3-2 excitation temperature, opacity and brightness temperature of N2D+ (dashed lines) and N2H+ (solid 
lines) on the left column, and of HCO"'" (solid lines), H^^CO+ (dotted line) and DCO+ (long dashed) on the right column. The column 
density is the same for all molecules and is 10^^ cm~^/(km pc~^). 
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Fig. 18. — Comparison of excitation temperatures (left) and line intensity ratios (rigth) of ortho-H2CO and para-D2CO for different 
transitions of both species and a common column density of 10^^ cm~^. Note that the quantum numbers of the selected transitions 
correspond to different symetries of the two isotopologues. Left panels show that excitation temperatures are different for both species for 
the range of densities found in dense cores. 



